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Abstract 

FASTERD is a Monte Carlo event generator to study the final state radiation 
both in the e^e~ 7r~^ir~^ and e~^e~ tt^tt^^ processes in the energy region 
of the 0-factory DA^^NE. Differential spectra that include both initial and final 
state radiation and the interference between them are produced. Three different 
mechanisms for the ttttj final state are considered: Bremsstrahlung process (both in 
the framework of sQED and Resonance Perturbation Theory), the (j) direct decay 
(e~^e~ (/o; /o + cr)7 — > Trir'j) and the double resonance mechanism (as 

e+e" p^vr^ — > 7r+7r~7 and e+e~ — > p — > un^ — > tt^tt^^). Additional models 

can be incorporated as well. 
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1 Introduction 



The anomalous magnetic moment of tlie muon (a^) is one of tlie most precise 
test of the Standard Model pQ. Theoretical predictions differ from the exper- 
imental result for more than 3a [2]. The main source of uncertainty in the 
theoretical prediction comes from the hadronic contribution, a^^"''^^ [2]. This 
contribution cannot be reliably calculated in the framework of perturbative 
QCD (pQCD), because low-energy region dominates, but it can be estimated 
by dispersion relation using the experimental cross sections of e+e~ annihila- 
tion into hadrons as an input [3] . About 70% of the hadronic part of the muon 
anomalous magnetic moment, a^^"''^\ comes from the energy region below 1 
GeV and, due to the presence of the p-meson, the main contribution to a^'*'^'^-* 
is related with the 71~^t:~ final state. 

Experimentally, the energy region from threshold to the collider beam en- 
ergy is explored at the ^-factory DA<I>NE, PEP - II and KERB at T(4S')- 
resonance using the method of radiative return (for a review see ^ and refer- 
ences therein). This method relies on the factorization of the radiative cross 
section into the product of the hadronic cross section times a radiation function 
H{q'^,9max,dmin) kuowu from Quantum Electrodynamics (QED) [5ll6|7|l8ll9] . 
For two pions final state it means that, in the presence of only the initial state 
radiation (by leptons, ISR), the radiative cross section a'^'^'^ corresponding to 
the process 

e+(p+) + e-(p_) ^ 7r+(pi) + n'ip^) + ^(k), (1) 



can be written as da'^'^"' = da'^'"{q'^)II{q'^,9max,9min) [HEE], where 9min and 
9max are the minimal and maximal azimuthal angles of the radiated photon, 
q = pi + P2 and the hadronic cross section a'"'^ is taken at a reduced CM 
energy. The final state (FS) radiation (FSR) is an irreducible background in 
radiative return measurements of the hadronic cross section [7p0] and spoils 
the factorization of the cross section. In any experimental setup the process of 
FSR cannot be excluded from the analysis. The KLOE experiment has devel- 
oped two different analysis strategies: the first one is with the photon emitted 
at small angle {6^ < 15°) and the other one is for the photon reconstructed at 
large angle (60° < 6^ < 120°), being for both 50° < 9^ < 130°. In the case of 
the small angle kinematics the FSR contribution can be safely neglected, while 
for the large angle analysis it becomes relevant (upto 40% of ISR). The large 
angle analysis allows to scan the pion form factor down to the threshold 

Radiative corrections (RC's) related to initial state radiation, i.e. the function 
H, can be safely computed in QED. For the FSR process the situation is dif- 
ferent. In the region below 2 GeV the pQCD is not applicable to describe FSR 
and calculation of the cross section relies on the low energy pion-photon inter- 
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Fig. 1. 

action model. Thus the measured FSR cross section gives an unique possibility 
to get very interesting information on the dynamics of interacting mesons and 
photons, to test the pion-photon interaction models and extract their param- 
eters pjj. 

In the case of neutral pions in the final states 

e+(p+) + e-(p_) ^ 7rO(pi) + 7r0(p2) +7(A;), (2) 



the ISR contribution is absent^], and the cross section is determined solely by 
the FSR mechanism. This process, together with asymmetries and cross 
section in the charged channel, allows to extract information on pion-photon 
interaction and test effective models for FSR. 

For realistic experimental cuts on the angle and energy of the final particles 
the cross section cannot be evaluated analytically and one has to use Monte 
Carlo (MC) event generator. The first MC describing the reaction ([T]) was EVA 
[7|. EVA simulates both ISR and FSR processes for non zero angle emitted 
photon [9^ > Omin)- For FSR the sQED model was chosen. Afterwards the MC 
PHOKHARA was written to include different charged final state and radiative 
corrections to ISR [14j. The contribution of the 0-meson direct decay, that 
is relevant at the DA$NE energy, was added, firstly in EVA [TJ] and then in 
PHOKHARA [13]. 

The computer code FASTERdS presented in this paper is a Monte Carlo event 
generator written in FORTRAN that simulates both processes ([1]) and ([2]), 
where the hard photon 7(/c) can be emitted by the lepton^ and/or the pions 
Fig. [H This program was inspired by EVA. The present version of FASTERD 
includes different mechanisms for the 71717 production: final BremsstrahlunJ^ 
in the framework of both Resonance Perturbation Theory and sQED, the con- 

^ This statement is valid only if one neglects multi-photon emission 
3 FinAL STatE Radiation at DA^>NE 

Only for the charged channel 
^ Only for the charged channel 
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tributions related to the intermediate state and the double vector resonance 
part (for details see Section [3]). Up to now only the case of one photon emis- 
sion has been considered. The code contains two types of source files: (1) the 
main program fasterd.f, where the calculations are done, and (2) the input 
file cards_fasterd.dat which defines the parameters for the generation. Both 
these files will be described in the following. As output of the program, the 
cross section of the process is evaluated and a PAW/HBOOK ntuple is pro- 
duced with the 4-momenta for each event of the outgoing particles. For the 
convenience of the user we have added a Makefile and the test output files for 
the Journal library. 

This paper is organized as follows. In Section [2] the main formulae for ISR 
and interference between ISR and FSR are presented. In Section [3] we give 
a general description of FSR process and present the FSR models that are 
included in our program. In Section H] the spectra for different FSR models 
are compared with analytical results and with MC PHOKHARA. Also a pos- 
sible generalization applicable to a wider region is described in Section |H In 
Section 0, we summarize the general procedure for calculating the spectrum. 
In Appendices A, B, C a short description of the input and output files is 
presented. 



2 Initial state radiation models and pion form factor 



The cross section of the processes ([T]) and ([2]) can be written as 



<2 



= C,,N\M\WdWdn- , \N = — / ^l - ^ j (3) 

where a is the fine structure constant, m.„ is the pion mass, u is the photon 
energy, Q = + s = and the invariant amplitude squared, averaged 
over initial lepton polarizations and summed over the photon polarization£3 
is 



|M|2 = \MisR? + \Mfsr? + 2Re{MisRM*psR)- (4) 



M^^^^^ (M*^^"^^^) corresponds to the ISR(FSR) production amplitude. The 
factor Ci2 = I for vr^vr" in the final state and Cu = 1 for ti'^ti'. 

6 We use J2poiar. e^ea = -9pa 
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For ISR process the invariant amplitude squared, averaged over initial lepton 
polarizations and summed over the photon polarizations, is 



|Ma^«)|2 = -A|F.(g2)|2^, (5) 



m 



q2 



ti \q' ) n )' tit2 

where rrie is the electron mass, xi,2 = 2p_,+ ■ p2, ti = —2p_k, t2 = —2p^k. 

The nonpoint-like behaviour of pions is determined by the form-factor (FF) 
F.jriq'^) that is the function of the pion mass squared g^. In the case of the 
neutral channel M^^^^^ = 0. Four different parameterizations for the pion FF 
are considered: Kiihn-Santamaria (KS) |l16J, Gounaris-Sakurai (GS) [17], the 
RPT parametrization and an "improved" version of Kiihn-Santamaria (see 
below) [H]. 

2.1 Kuhn-Santamaria and Gounaris-Sakurai pion FF 



Based on the results of Ref. [T6j, the pion FF describing the p — uj mixing and 
the first excited p resonance (p'), can be written as 

F.(,^) = + (6) 



where for the KS parametrization 

2 

ml — — i\f(pTj.{q^) 
and for the GS one [17] 



ml + if(0) 



ml — + H{q^) — i^/q^Tr{q'^ 



with 



H{q' 



pliml) 
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The energy dependence for the p mesons is taken in the form 

rp(g^) = rp^f^^V-0(^^-4m^)- (9) 



For the uo resonance a simple Breit-Wigner resonance form with constant width 
was used for both parameterizations. 

As one can see the single resonance contribution is normalized to unity at 
q^ = for both parameterizations (-Br(O) = 1) whereas the right normalization 
for the pion FF, -F^(O) = 1, is realized by the corresponding choice of the 
parameters a, [3. 

All model parameters (the mass and width of the resonances as well as the 
parameter a, (5) are determined in the input file cards_fasterd.dat. The KS 
pion FF parametrization corresponds to the function F_pi_ks whereas the GS 
one to F_pi_gs. 



2.2 RPT parametrization for the pion FF 



The Resonance Perturbation Theory is based on Chiral Perturbation The- 
ory (xPT) with the explicit inclusion of the vector and axial-vector mesons, 
Po(770) and ai(1260). Whereas xPT gives correct predictions for the pion form 
factor at very low energy, RPT is the appropriate framework to describe the 
pion form factor at intermediate energies {E ~ nip) [12]. According to the 
RPT model the pion form factor, that describes the p — uj mixing, can be 
written as: 

FM') = 1 + ^4^f'(^')fl - ^B^^'i^'A (10) 
In \ 3m^ J 



where q is the virtuality of the photon, = 92.4 MeV and the parameter 11 
describes the p-uj mixing. As before a constant width is used for the cu-meson. 
We also assume that the parameter IIp^ is a constant and is related to the 
branching fraction Br{uj tt+tt"): 

5r(c. ^ TT+vr-) = -ffl^. (11) 
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As was mentioned above, the value of Fy and Gy, as well as the mass of 
the p and u mesons {nip and m^, correspondingly), the parameter of the p-u 
mixing IIpj^ and the width of the uj meson are determined in the input file 
cards_fasterd.dat. The RPT parametrization corresponds to the function 
F_pi_rpt. 

Inclusion of the p' meson modifies the form factor (fTOl) as 



FAf) = l + ^^B-Hf) (l - (12) 



2.3 "Improved" Kuhn-Santamaria 



In Ref. the pion FF was estimated in the framework of the dual-QCDTv^^oo 
model. However this model assumes the resonances to be of the zero width. 
The authors of Ref. |i8| included the final width of the p° and the three lowest 
excited p' states and obtained the following expression for the pion FF 

FM') = E c„i?f ^(g^) + E cn^^. (13) 

n=0 n>4 ^ 



This parametrization for the pion FF is contained in the function F_pi_ks_new. 
The numerical value for the parameters c„ and m„ (the mass of the excited p 
mesons) is taken from Ref. [I8] . 



3 Final state radiation models 



Using the underlying symmetry, like gauge invariance, charge-conjugation 
symmetry of the final particles and the photon crossing symmetry, it is pos- 
sible to write the FS tensor M^'^\ that describes the 7* — >■ 7r+7r~7 vertex, in 
terms of three gauge invariant tensors (see [21] and Ref. [23,24] therein): 



Mr(Q, k, I) = rr/i + rr/2 + r^fs, (14) 

Tr = k^Q''-g'"'k-Q, I = pi-P2, 

T^" = k ■ lil^'Q'' - g'^'k ■ I) + rikf'k -l-l^k- Q), 

T^'' = s{g^"'k ■ I - WV) + Q^{rk -Q-Q^k- 1). 
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The model dependence comes in only via the implicit form of the scalar func- 
tions fi (we will call them structure functions). 

Thus the FSR and interference {ISR * FSR) part to the invariant amplitude 
squared (jlj) is 



M, 



FSR\ 



aii|/ir + 2ai2Re(/i/*) + a22|/2 



2a23i?e(/2/3*) + agsl/sr + 2anRe{hf^) 



(15) 



and 



Re{MisRMpsR) 



(16) 



+ A,ReiF^iq')f;] 



The form for the coefficients and Ai can be found in Ref. [21] Eqs. (17), 
(26), correspondingly. 



Here is the list of the FSR production mechanisms that are included in FAS- 
TERD: 



TT^ + TT + 7 Bremsstrahlung process (17) 

~^ (/o; /o + cr)7 — TT + TT + 7 (f) direct decay (18) 

(0; uj') -^p7r-^7r + 7r + 7 Double resonance process (19) 
(p/p') — > ton^ 7r° + 7r° + 7 Double resonance process (20) 

Thus the total contribution to the functions fi introduced in f[T^ is 

/. = f^'""^ + /i'^ + ft"'^ (21) 



e+ + 
e+ + e" 



In the next section we present the models describing these processes. 



3. 1 Final state Bremsstrahlung 



Usually the combined sQED*VMD model is assumed for the FS Bremsstrahlung 
process [7111^ . In this case the pions are treated as point-like particles (the 
sQED model) and the total FSR amplitude is multiplied by the pion form fac- 
tor Fff(s), that is estimated in the VMD model. Unfortunately the sQED*VMD 
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model is an approximation that is valid for relatively soft photons and it can 
fail for high energy photons, i.e near the vr+vr" threshold. In this energy region 
the contributions to FSR, beyond the sQED*VMD model, can be important. 
As mentioned in Section [21 RPT is supposed to be an appropriate model to 
describe the pion-photon interaction in the region about and below 1 GeV and 
this model is used to estimate the contributions beyond sQED*VMD. 

Using the sQED*VMD model the structure functions fj-^"^^"^^ (^^ee Eq. (l?I|) ) 
are 



sQED _ 



fl 

fsQED 



(k-l) 



2' 



sQED 



-2FJs) 



{k-Qf-{k-l) 



2 ' 



(22) 
(23) 



In the framework of RPT the result is 



where 



n 



Mr 



Fy — 2FvGy 



2 + 



D{l)D{-l) 



1 

— + 

mi 



m' 



- s - iv/ir^ 



+ 



[s + k-Q)[Aml-is + l^ + 2k-Q)] 



Ami 



8D{l)D{-l) 
s + l^ + 2k-Q) 



flml 8D{l)D{-l) 



k-l 



f^ml2D{l)D{-l) 



D{1) 



ml-{s + f 



(25) 
(26) 

2kPQ + m)/A.{27) 



For notations and details of the calculation we refer the reader to [21]. The 
functions A//^^"^ are calculated by the subroutine chpt_rpt whereas //*^^^ 
are evaluated in the functions fsrl_sqed, fsr2_sqed, fsr3_sqed. 

We would like to mention here that the contribution of any model describing 
Bremsstrahlung FS process can be conveniently rewritten as in Eq. (!24l) and 
in the soft photon limit the results should coincide with the prediction of the 
sQED*VMD model. 
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3.2 direct decay 



As it was mentioned in the Introduction, at the DA$NE energy (a/s = 
= 1.01944 GeV) the final state Trn'-f can be produced via the intermediate 
meson state. In this section we consider the direct rare decay — >■ irn'j. As 
shown in [T5|22] this process affects only the form factor fi of Eq. ([1] 

ml + im^T^ ' 



The direct decay is assumed to proceed through the intermediate scalar 
meson state: — (/o + o")7 — nnj and its mechanism is described by a single 
form factor 

In the input file the different following models can be chosen to describe this 
decay: 



• Non structure model 

• Linear sigma model 

• Chiral unitary approach [25] 

• Achasov kaon loop model 

• Achasov kaon loop model with inclusion of the a meson [27| 

• Non structure model that includes both the /o and a mesons [22] 

The explicit expression for the factor /^(s) for the first three models can be 
found in p^. For the fourth ("Achasov kaon loop") model the form factor 
reads: 



2Ti'^m\\ 


>\ 


- s + Relif^Xm), 


)-%(Q^)) 




s 







(29) 



where /(., .) is an analytical function [T5|28] and the phase 6b{s) = 6y s — ^m^, 
b = 75°/GeV [27]. The term ReIlfQ{m'j^) —11^2(5) takes into account the finite 
width corrections to the /o propagator [26] . 

In a refined version of this model which includes the a meson in the interme- 
diate state [27], the form factor can be written as 



j^K+K- (^^^-^ ^ g<t>K+K- ^i(5^,(s)+SKK(Q'')) , 



m 



2 



11 



R,R' 

where Grr' is the matrix of inverse propagators [27] . Such an extension of the 
model improves the description of the data at low mTrvr and was recently used 
by KLOE in the fit of ^ n^n^'j spectrum [29] . 

The contribution of the direct decay is calculated in the subroutines fphl_bcg, 
fphl_In, fphl_cpt, fphl_a4q, fphl_a4qs, fphl_pac according to the cor- 
responding model. 



3.3 Double resonance contribution 



Another mechanism producing the final inT'-f state is reported in Eqs. (fT9|) . 

(EQD. 

First we consider the (j) pn mechanism. In this case the meson decays 
in [p^ir'^) for the charged channel and in [p^T^^) for the neutral one and then 
p — >• 717. The corresponding contribution to the functions Z/'"^'^*^ of Eq. (!2T|) is 



ft' = Y^l(^-Q + l')iDp{Rl) + D,{Rl)) + 2k ■ l{D,{Rl) - D,{Rl))l 
fr = ^[D,{Rl)~D,{Rl)l 



where 



m% - s - im^T^ 



and Dp{R) = ml - R - impy/RTp{R), = {k + pi)^, R_ = {k + p2?- The 
quantities g'^^, gl^^ are the coupling constants determining respectively the 



pvr and p —>■ tt'j vertexes correspondingly, 



/3r(0 e+e- 



To make connection with the KLOE analysis [29] we add also the phase of the 
uj-cj) meson mixing and the constant factor nj^*^^ 



^ Including H^*^^, in our opinion, rescales the constant gfn that cannot be directly 
determined from any experimental decay width 
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Thus the coefficient (1301) is rewritten as 



CM - ^S^nr-e-.. (31, 



In the energy region of DA<I>NE also the tail of the excited u meson can play 
a role: 7* — cj' — > pvr. The explicit form of this contribution is written similar 
to (EOD and is: 



VA^ gp^ g'^l F^. 



s 



In the considered energy region Eq. (]32l) can be approximated by a complex 
constant C^^ whose the numerical value is taken from the fit of tt^tt^'j spec- 
trum 1,291 



Therefore we have 



^(»)-S^nr-e-- + C^. (33, 



For the neutral case the 7* {p/p') ^^^^ 7r°7r°7 reaction contributes as 
well. As for the uj' pvr channel we replace the explicit expression 

_ 9L Fp' s VA^ g':^^ s 

'^ujn — 2 '■ T=i ' 9 '■ \'^^) 

mp — s — inipil pi nip — s — tmpL p 



by a complex constant C^,^, which is again taken from the fit. Then the contri- 
bution of the 7* —>■ {p/p') —>■ onT^^ 7r''7r°7 mechanism to the functions f^^^'^^^ 
is 



CP 

fr = TT^iik ■ Q + l^){D^{Rl) + D^{Rl)) + 2k ■ l{D^{Rl) - D^{Rl))l 

iOTTO; 



CP 
lovra 



CP 

f3'^=Y^mRi)-D^{Ri)]. 



Also if we suppose that the direct decay 7* — > /ovr ^ tttt'j takes place it contributes 



to C;,. 
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Finally the total contribution is 

^i.ect)^^iM^.S,^^fip^)^ (35) 

where c = 1 for the neutral final state and c = for the charged one. We 
include also an additional phase between the double resonance and (j) direct 
contributions (5p). 

The subroutine rhotopig calculates the contribution for the double resonance 
mechanism. 



4 Comparison of MC results with analytical prediction. Discussion 
of the results 



For the full angular range 0° < 9^, Oj, < 180° it is possible to obtain the 
analytical expression of the cross section for the processes ([T]) and ([2]). Here 
present the results for the charged final stat^. 



we 



The ISR cross section is 



whereas for the FSR process we have 



(2/ii-s/i2), (37) 



where the form of the functions hi and the equations to calculate them using 
the functions fi can be found in Ref. |2lj. For example, in the case of sQED 



167re|F.(.)P 



ikQr+r--ml\i-q' + iq'-2ml)L,) 



{31 



/ 4m^ 1 1 + ^ 

with ^ = \ I TT and Li = - In -. 

V r ^ l-'^ 



^ For the 7r''-7r'^7 case the contribution both from the (p direct decay and the double 
resonance vector mechanism is a half of the corresponding one for the 7r"'"7r~7 final 
state (except for the 7* ^ p ^ lott case that appears only for the 7r'^7r'^7 FS). 
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\ (GeV^) 



Fig. 2. Comparison of the analytical and numerical results for ISR: 0' < < 180', 
0' < 0^ < 180'. Left: the solid line corresponds to the analytical result whereas the 
points with errors have been obtained by FASTERD. Right: the ratio between the 
analytical results and numerical ones is fitted by the constant ^o- 



Fig. 3. Comparison of the analytical and numerical results for different contributions 
to FSR: 0' < 6'7r < 180', Q' < 6^ < 180'. For the (f) direct decay the kaon loop model 
with /o + fj is chosen whereas the RPT parametrization is chosen for the pion FF. 
The solid lines correspond to the analytical result whereas the points (triangles, 
squares) have been obtained by FASTERD. 

In Figs. [2] and [3] we compare the analytical predictions fl36|) and fl37j) with the 
numerical results obtained by FASTERD for the cross section of the reaction 
e^e~ '7r+7r~7. The ratio between the analytical and numerical results is 
fitted by the constant Aq. We have good and the value of Aq compatible 
with one. 

An important feature of FASTERD is a possibility to estimate the double 
resonance contribution. At the first sight in the energy region about s ~ 
it is enough to include only the 7* ^ — >■ pvr — tttt'^ mechanism. However, 
at s = m\ our simulation gives the ai^cp — » ((/o + cr)7 + p°7r°) — > 7r°7r°7) = 




15 




q" (GeV") q- (GeV^) 



Fig. 4. Different contributions to FSR for the 7r'^7r''7 (left) and 7r+7r^7 (right) final 
states: 0' < 6*^ < 180' , 0' < 9-y < 180' . The filled circles are for the double resonance 
contribution, the triangles correspond to the (/o + (7)7 mechanism and the empty 
circles correspond to the interference term between the (/o + cr)7 and the double 
resonance contributions. 



0.451 ± 0.001 nb whereas cj((p/p') ^ un^ 7r°7r°7) = 0.529 ± 0.012 nb 
(in agreement with KLOE results [30]). In the case of the 7r+'7r~7 the double 
resonance mechanism (only the — >■ p^vr^ contributes) does not give a large 
contribution in the high region whereas it gives visible the cross section in 
the low region (Fig. HI right). In the current version FASTERD allows to 
simulate the processes ([T]) and ([2]) for s ~ m|. To simulate the processes for 
lower (upto the threshold) and higher (say upto 1.2 GeV^) energies the exact 
value of the constants C^^ and C^^ (Eq. fl32|) and Eq. flMl) ) has to be used. 

We compare the ISR spectrum with the results from PHOKHARA [14J. The 
'improved' Kiihn-Santamaria parametrization was chosen for the pion form 
factor and PHOKHARA was run at the leading order approximation (only 
one photon is radiated). The results of comparison are presented in Fig. [51 
There the ratio between FASTERD and PHOKHARA results is fitted by 
the constant Aq. As one can see the PHOKHARA result coincides with the 
FASTERD prediction. As it was discussed in Ref. [H] the radiative corrections 
due to an additional photon emitted by the leptons are relevant and we are 
planning to include them in the code. An additional photon can be emitted 
either from initial state or from final state. Inclusion of the second photon 
emitted from the initial state radiation is a technical task and can be obtained 
by replacing the leptonic tensor (it has been calculated already [THI2T] ) and 
substitution k ^ k' , Q ^ Q' = — — k' in Eq. (fl^ . where k' is the energy 
of the additional photon emitted from the initial state. Also the one photon 
phase space has to be replaced to the two photon phase space. For the second 
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I- 87.22 / 98 
Ad 0.9995 ± 0.1392E-02 



il i l iii l i li ii jit ii ttt i it , # i VitJitll i ti i llti i i i t i i'i^^ 



(GeV") 

Fig. 5. Comparison of the numerical results for ISR simulated by FASTERD and 
PHOKHARA :0 <9„< ISO', 0' < 9^ < ISO'. The "improved" Kiihn-Santamaria 
parametrization is chosen for the pion FF. The ratio between the numerical results 
obtained by FASTERD and PHOKHARA is fitted by the constant Aq. 

photon from the final state the situation is more difficult as there is not yet a 
model-independent calculation for the final state tensors with two photons. 



5 The calculation of the spectrum. General algorithm 

The logical structure of FASTERD is presented in Fig. El 

Input parameters are read from the file cards_fasterd.dat that is modified by 
the user (see Appendices A, B for the description of the input file) . Here it is 
possible to choose the name of the output files. The main part of the program 
runs twice. On the first step, in order to estimate the maximum value of the 
integrand (Mmax) some trial events (NMPT) are generated. On the second 
step (NEVE) events are produced and the acceptance-rejection method is 
applied to evaluate the cross section da/dq^. The generator RANMAR [3T] is 
used to generate the random numbers. The variables (g^, cos6'-y, cos0^) are 
generated in the laboratory frame {e^e~ center of mass system), the variables 
(cos^^ and coscf)'^) are generated in the pion mass reference system. They are 
calculated in the subroutines qquadrat and photonangles. The subroutine 
vectorb transforms the azimuthal and polar angles of pions in the laboratory 
system. The subroutine rejection checks the restrictions imposed by the input 
file. The matrix squared element is produced inside the function mat. Using 
the value mat returned by this function the differential cross section inte is 
computed as inte = mat*num* jacjpion* jacjphot, where num is a kinematic 
factor for the cross section, jacjpion and jacjphot are the replacement factors 
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(they are calculated in the subroutines qquadrat and photonangles). The 
information about the accepted events is stored in the file output. hbk. At the 
end of the generation, the cross section is computed by the ratio of accepted 
to generated events, and the error is estimated. The information about the 
number of generated and accepted event as well as the value of the cross 
section and its error is written in a output file. 

Function mat. The input for mat is the 4-momenta of all particles in the 
laboratory frame and the type of process (ISR, FSR or both ISR+FSR). Ac- 
cording to the type of the process either the subroutine matrJsr or matr_fsr 
or matrjnt are called. They compute the Eqs. (IT^ . (IT^ presented in 
this paper. 

Subroutine qquadrat. The following probe function is used to generate the 
variables g^, cos^tt and coscpT, (we remind that the angles of the pions are 
generated in the pion reference system) 



' s{s — q^) ~^ (g^ — m?pY + F^m^ ~'~ 

1(1 1 



s^s-g-j Vl- ^1- ^cos^^ l + ^l-^cos^^y 

The polar pion angle is generated uniformly. 

It results in the following expression for the jac-pion 

2nV , , 

iacjpion = — — , (39) 

/(g^cos6'^) 



where 



V = — m H arctan — arctan - 



:ln 



^ , /I 4m2 



sil - ^/l - ^)V1 - ^ 1 - VI - ^ 
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Input: choice of: 

- final state ('7r'^7r°7 or 7r+7r^7); 

- type of process (ISR (7r+7r-7), FSR, ISR+FSR); 

- experimental parameters and cuts: s, 6.,^.^^{min;max) etc.; 

- pion FF (7r+7r^7), FSR model and parameters 



First run: the max value 
of integrand M^ax 




Generation of 4-momenta: 

Kinematics selection check: 

6T,{min) < 9^ < 9^(max) etc. 
Matrix squared element calculation 




t 


Second run: the cross section 
da /dq^ 







Output files: 
output. hbk: event information 
output. out: da /dq^ value etc. 




For the second run: 

- storage of the accepted events 

- calculation of the cross section 

da/dq^ and uncertainty 





Fig. 6. The logical structure of FASTERD 
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Appendix A. Input file 



The main parameters are set in cards_fasterd.dat input file. It defines the 
cuts used in the generation, and the parameters needed for the for the FSR 
model and the pion FF. 
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The following variables are included in cards_fasterd.dat file: 

• NEVE - number of the generated events to be generated 

• NPTM - number of events to find the maximum of the cross section. 

• NRND - random seed 

• FSKIND - kind of final state (FSKIND=0 for 7r°7rO and FSKIND=1 for 

• KIND - two parameters which allows to choose the kind of radiation: (ISR 
correspond to (1 0), FSR to (0 1) and ISR+FSR to (1 1)) 

• FPKIND - defines the pion FF parametrization (1 for KS, 2 for GS, 3 for 
the pion FF in the framework of RPT and 4 for "improved" KS) 

• RHOKIND - choice of FSR contribution 

• BREMKIND - choice of Bremsstrahlung model 

• fOKIND - choice of the model for the (f) direct decay 

• RESPIKIND - allows to choose the intermediate states (pvr and/or uj-k^) for 
the double resonance contributiorF°1. 

• HOUT - defines the name of the output file output. hbk where information 
of the generated particles (4-momenta) is inserted. 

The parameters that define the experimental cuts are 

• ENES - the squared energy of the initial particles 

• GMIN - the minimal energy of the radiated photon 

• GMAX - the maximal energy of the radiated photon 

• ACUT - the minimal and maximal azimuthal angles of photon and pions 

The value of the constants for the pion FF are defined in the arrays FP 
and FPIHIGH. The remaining set of arrays FOpar,FOfase,FOnonstr, VMDPAR 
defines the numerical values for the direct decay and the double resonance 
contribution. 



Appendix B. Example of input file 

As an example, the input file cards_fasterd.dat is set to generate the e'^e 
7T+7c~'y process where: 

• the photon is emitted both from IS and FS; 

• the RPT parametrization is chosen for the FF; 

• FSR is given by sQED Bremsstrahlung and (f) —>■ {fo + a) decay; 

• the experimental cuts are: 50° < ^7r,7 < 130° 

• the output file is named isr_rpt_fsr_sqed_f0+sig_50_130_2e6.hbk 



The UJTT intermediate state is produced only for the final state 
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C* number of events to be processed 

NEVE 1000000 
C* — random seed 
NRND 9513463 

C* number of points used to evaluate the maxima 

(at least 10**4) 

NPTM 10000 

C* — piOpiO (FSKIND=0) or pi+pi- (FSKIND=1) channel 
FSKIND 1 

C* type of radiation 

= (1 D) 



initial state radiation on 
final state radiation on 
initial+f inal state radiation 



C*KIND 1 
C*KIND 1 
C*KIND 1 1 
KIND 1 1 

C* Model for the Pion form factor 

KS 
GS 
RPT 

KS+dual-QCD 



1 
2 
3 
4 



C*FPKIND 
C*FPKIND 
C*FPKIND 
C*FPKIND 
FPKIND 3 

C*FSR brems 

rhokind 1 
C*bremkind 1 
C*bremkind 2 
bremkind 1 
C*fOKIND 1 
C*fOKIND 2 
C*fOKIND 3 
C*fOKIND 4 
C*fOKIND 5 
fOKIND 5 
C*respikind 1 
C*respikind 2 
C*respikind 3 
respikind 3 

C* Parameters of the Pion Form factor 

C*RPT MRHO GAMMARHO MRHOL GRHOL MOMEGA GOMEGA AL BE argA 
pi\_rho FV FVl GVl 

FP 0.774 0.143 1.37 0.51 0.7827 8.68E-3 0. 0. 0. 
-0.0027 0.154 0. 0. 

C*FPIHIGH c_0_pion c_2_pion c_3_pion m_rho2_pion g_rho2_pion 
m_rho3_pion g_rho3_pion 

FPIHIGH 1.171 0.011519 -0.0437612 1.7 0.24 2.0517 0.41 

C* Parameters of Scalar contribution (f 0,f 0+sigma) in FSR 



fO phirhopi 
1 
! sQED 
! RPT 

linear sigma model 
non structure model 
chiral unitary approach 
kaon loop with fO only 
kaon loop with fO+\sigma 

only \rho\pi\gamma 
only \omega\pi\gamma 

both \rho\pi\gamma and \omega\pi\gamma 
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C*FO+SIGMA GfO\_k+k- GfO\_p+p- GPHI phase (deg) fOMASS msigma 

gsigpp gsigkk CfOsig 

FOPAR 4.02 -1.76 4.482 0. 0.982 415.0e-3 -2.2 -0.37 0.015 
C*fOphase FO+SIGMA mOk m2k LambdaK bOp blp b2p Lambdap 
FOfase 0.363 0.757 1.24 5.4 3.7 5.0 0.200 
C*F0nonstr gfO_pp gf0_kk gphi_fOg gphi_sgg aO 
FOnonstr 0.57 1.14 3.5 5.0e-3 1. 

C* Parameters of double resonance mechanism 

C*VMDPAR g\_rhopig g\_phirpi prhores beta\_bro beta\_wphi 

c\_modrpi c\_pharpi c\_modwpi c\_phawpi 

VMDPAR 0.295d0 0.811d0 0.677d0 32.996d0 163. OdO 

0.26d0 3.1112d0 0.85019d0 0.45d0 

C* histo output file 

*HOUT 'isr_rpt_fsr_sqed_f0+sig_50_130_le6.hbk' 
ENES 1.039202865 
GMIN 0.02 
GMAX . 600 

C*ACUT theta_min_\gamma theta_max_\gamma theta_min_\pi theta_max, 
ACUT 50. 130. 50. 130 

C* 

STOP 
END 



Appendix C. Output file 

There are two kinds of output files: (output. hbk, HBOOK/PAW ntuple for- 
mat, defined in cards_fasterd.dat) contains the 4-momenta of outgoing par- 
ticles for each accepted event; (output. out, text format, output of the pro- 
gram) contains the number of generated and accepted events, the value of the 
cross section with the error, and the biggest value of the integrand. 
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